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The structures of embedded clusters in the Perseus. 
Serpens and Ophiuchus molecular clouds 
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ABSTRACT 

The young stellar population data of the Perseus, Ophiuchus and Serpens molecu- 
lar clouds are obtained from the Spitzer c2d legacy survey in order to investigate 
the spatial structure of embedded clusters using the nearest neighbour and minimum 
spanning tree method. We identify the embedded clusters in these clouds as density 
enhancements and analyse the clustering parameter Q with respect to source lumi- 
nosity and evolutionary stage. This analysis shows that the older Class 2/3 objects 
are more centrally condensed than the younger Class 0/1 protostars, indicating that 
clusters evolve from an initial hierarchical configuration to a centrally condensed one. 
Only IC 348 and the Serpens core, the older clusters in the sample, shows signs of 
mass segregation (indicated by the dependence of Q on the source magnitude) , point- 
ing to a significant effect of dynamical interactions after a few Myr. The structure of 
a cluster may also be linked to the turbulent energy in the natal cloud as the most 
centrally condensed cluster is found in the cloud with the lowest Mach number and 
vice versa. In general these results agree well with theoretical scenarios of star cluster 
formation by gravoturbulent fragmentation. 

Key words: stars: formation - ISM: clouds - ISM: kinematics and dynamics - open 
clusters and associations: general - infrared: stars - methods: statistical 



1 INTRODUCTION 

It is believed that stars form from high-density regions 
in turbulent molecular clouds that become gravitation- 
ally unstable, fragment and collapse. That way, a star 
cluster is built up in a comple x interplay between grav - 



ity and supersonic turbulence ([Mac Low & Klessen 20041 : 
iBallesteros-Paredes et all l2007l : McKce & Ostrike^ 
Therefore the structure of an embedded cluster, i.e. the 
spatial distribution of its members, may hold important 
clues about the formation mechanism and initial condi- 
tions. Several characteristics of star-forming clusters seem 
to be correlated with the turbulent flow velocity, e.g. 



lar medium shows a hierarchical structure (sometimes de- 
scribed as fractal) from the largest GMC scales down to 
individual cores and clusters, which are sometimes hierar- 
chical themselves. There is no obvious change in morphology 
at the cluster boundaries, so clusters can be seen as the bot- 
tom par ts of the hierarchy, where stars have had the chance 
to mix (lEfremoy fc Elmegreen|[l99i : lElmegreen et al.ll2000l : 



J, ro mix (li^iremoy 

lElmegreenll2006l ). 



the star formation efficiency llKlessen. Heitsch fc Mac Low 
200C ; IClark fc Bonnellll20o4 ISchmeia. Klessen fc Froebrich 
200E ) , mass accretion rate s ISchmeia fc KlessenI |20 04|) , or 



the core mass distribution (|Ballesteros-Paredes et al."2006l) 
and the initial mass function (Padoan fc Nordlund 2005). 

Giant molecular cloud (GMC) complexes usually con- 
tain multiple embedded clusters as well as more distributed, 
rather isolated young stellar objects (YSOs). The interstel- 
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Since embedded stars are best visible in the infrared 
part of the spectrum, and since nearby GMCs can span 
several degrees on the sky, wide-field near-infrared (NIR) 
and mid-infrared (MIR) surveys are useful tools for identi- 
fying and mapping the distribution of young stars in GMCs. 
The advent of the Spitzer Space Telescope has in particular 
led to large advances in such studies. The S pitzer Cores to 
Disks (c2d) Legacy programme survey ([Evans et al.ll2003l ) 
aimed at mapping five large nearby star-forming regions, 
including Perseus, Serpens and Ophiuchus using Spitzer 's 
IRAC (3.6 to 8 Mm) and MIPS (24 to 160 ^in) instr u- 
ments |j0r"gensen et all l2006l: iHarvev et all |2006| . l2007al lbl: 
iRebuU et al.ll2007l : IPadgett et al.ll20oi r 

The Perseus molecular cloud complex is an extended 
region of low-mass star formation containing the well- 
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studied embedded clusters IC 348 and NGC 1333. The 
latter is considered to be significantly young er (< 1 Myr; 
iLada. Alves fc Lad a"l996': Wilking ct al. 2004) than IC 348 
3 Myr: iLuhman et al.. .2003 : Mucnch ct al., .20071 ). The 
distance to the Perseus cloud is uncertain and may lie 
somewhere between 200 and 320 pc, presumably there is 
a dist ance gradient alon g the cloud (see e.g. the discus- 
sion in lEnoch et"aLll2006l ). Areas of 3.86 deg^ and 10.6 deg^ 
have been mapped by IRAC and M I PS, respectively 



Table 1. Classification scheme 



lave been mapped by IKAU and M i r a, respectively 
J0rgensen et al'lbood : ikebull et al]|2007^ . |j0rgensen et al.1 



20061) identified a total number of 400 YSOs, thereof 158 in 
IC 348 and 98 in NGC 1333, which means that a significant 
frac tion of young stars is fo und outside th e main clusters (see 
also iHatchell et al.|[20oH ). iMuench et al.1 (^007 ). also based 
on Spitzer MIR data, report a known membership of IC 348 
of 363 sources and estimate a total population of more than 
400. 

The Serpens molecular cloud is a very active star- 
forming region covering mor e than lOdeg^ on the sky at 
a dis tance of about 260 pc l|Straizys. Cernis fc Bartasiiitel 
1 19961 ). Areas of 0.89 deg^ and 1.5 deg^ have been mapped 
by IR AC and MIPS, respectively (|Harvev et al.1 |2006| . 
2007al). In th is regio n, 235 YSOs have been identified by 



Harvey et all (|2007bl ). The age of the main cluste r, called the 
Serpe ns cloud core, is estimated to be ~ 2 Myr l|Kaas et al.l 
|2004| ). 

The Ophiuchus (or p Ophiuchi) molecular cloud is one 
of the closest star-form ing regions at a distance of about 
135 pc iMamaiekl l2008h . Its main cluster, L1688, is the 
richest known nearby embedded cl uster. It is very young , 
with an age of probab ly < 1 Myr (jCreene fc Meven 1 19951 : 
iLuhman fc Riekelll999l ). 14.4 deg^ have been observed with 
MIPS, leading to th e identification of 323 YSO candidates 
l|Padgett et al.ll2008l ). 

In this paper we analyse the structure and distribution 
of the young stellar population unveiled by the c2d pro- 
gramme in these nearby molecular clouds and compare the 
results to current theoretical scenarios of star formation. For 
this purpose, we retrieved the c2d point source catalogues of 
the Perseus, Serpens and Ophiuchus molecular clouds and 
applied two complementary statistical methods to identify 
embedded clusters and analyse their structures. The selec- 
tion of the data set is described in Section [5] the statistical 
methods and their application to the data are explained in 
Section [3] and the results and discussion are presented in 
Section m and [S] 



2 DATA SELECTION 

The third delivery of data from the c2d legacy project 
lEvans et al.ll2005l ). which combines observations from the 
IRAC and MIPS cameras was obtained from the c2d web- 
sitfl The combined photometric catalogues for the Perseus, 
Ophiuchus and Serpens molecular clouds were retrieved. 
These catalogues are products resulting from extensive anal- 
ysis of the photometric observations in the four IRAC bands 
(central wavelengths 3.6, 4.5, 5.8 and 8.0 /im) and the three 
MIPS bands covering the 24 to 160 fim wavelength range. 



^ |http : //ssc . spitzer . caltech. edu/legacy/c2dhlstory .htmll 



object type 
CEvans et al. 2005) 



classification, 
main criteria 



red 

redl 

red2 

YSOc_MPl.red 
YSOcJR4-|-MPl.rcd 



Class 0/1 

(objects only visible at 
wavelengths ^ 8 iim) 



YSOcJR4 

YSOcMPl 

YSOcJR4-|-MPl 

YSOcJR4.PAH-em 

YSOcJR4_star+dust(BAND) 

YSOcJR4.red 

YSOc_MPl.PAH-em 

YSOc_MPl.star+dust(BAND) 

YSOcJR4-|-MPl.PAH-em 

YSOcJR4-(-MPl_star-|-dust(BAND) 



Class 2/3 

(objects visible at 
wavelengths < 8 /im, 
SEDs with slopes 
typical for Class 2/3) 



star 

star-l-dust(BAND) 



stars 

(SEDs fitted with 
photospheres) 



The sensitivity of the survey is expected to sample objects 
down to 0.1 - 0.5 M© in all three clouds. 

The delivery catalogues provide source identifications 
and classifications based on a consistent method of anal- 
ysis for all the targets covered by the legacy survey as 
described by Evans ct al. (2005). The catalogues contain 
positions, photometry and associated quality flags, source 
type such as point /extended source and source classification 
such as stars/galaxy /YSO. These catalogues focus on iden- 
tifying and classifying YSOs and distinguishing them from 
stars and galaxies. Extensive analysis of the point sources 
using colour-colour and colour-magnitude diagrams, spec- 
tral indices in the observed bands, comparison with SWIRE 
{Spitzer Wide-area Infrared Extragalactic survey) data to 
remove galaxy and PAH emission contaminants is used for 
this purpose. 

In the present work we use these catalogues to retrieve 
the stellar and YSO population in the three molecular clouds 
based on the source classification flags. Although YSOs are 
classified based on their spectral energy distributions (SEDs) 
and MIR colours, a certain amount of confusion is always 
present when the complete SEDs are not available and only 
a part of the SED is sampled. The IRAC and MIPS obser- 
vations cover the wavelength range where the YSOs lend 
themselves for an effective classification, however, confusion 
arising due to partial sampling of the SED cannot be re- 
moved. Further, an edge-on Class 2 source can imitate a 
Class 1 source etc. when classification is made using partial 
SEDs. Based on the nature of the point source SED, which 
is represented by 19 different flags in the c2d catalogues, 
we group the point sources into three categories, namely 
stars (foreground/background and cloud members), YSOs 
in Class 0/1 phase and YSOs in Class 2/3 phase. While 
the c2d catalogue provides a distinct classification for stars, 
the YSOs are represented by 15 different object types. We 
use the thumb rule of grouping sources that are visible only 
in wavelengths longer than 8/xm into the Class 0/1 group 
and the remaining into the Class 2/3 group. This leads to 
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the classificatio n given in Ta b le [J where we relate the ob- 
ject types from lEvans et alj (|2005l ) and their characteris- 
tics to the common evolutionary classes. Our classification 
does not necessary represent the 'true' nature of the objects, 
which would in some cases require data from wavelengths 
longer than provided by Spitzer and a more sophisticated 
SED modelhng. 



3 STATISTICAL METHODS 
3.1 Nearest neighbour method 

While clusters are readily identifiable as peaks in stellar den- 
sity maps, the exact delimitation of a cluster is difficult 
and will always be somewhat arbitrary, since all the em- 
bedded clusters are surrounded by an extended population 
of YSOs distributed throughout the entire clouds. Several 
methods to detect clusters in a field have been developed. 
Usually they are based on stellar densities (derived e.g. from 
star counts or Nyquist binning) and consider as clusters 
all regions above a ce rtain deviation (~ 2 — 4g) from the 
background level (e.g. Lada et al. 1991 : Lada fc Ladal 19951: 



Kumar. Keto fc ClerkinI l2006l : iFroebrich. Scholz fc Raftervl 
20071 ). 



The nearest neighbour (NN) method (|von HoerneJ 
ll963l : [c"asertano~fc Hut. 1985. ) estimates the local source den- 
sity pj by measuring the distance from each object to its jth 
nearest neighbour, where the value of j is chosen depend- 
ing on the sample size. There is a connection between the 
chosen j value and the sensitivity to those density fluctua- 
tions being mapped. A small j value increases the locality 
of the density measurements at the same time as increasing 
sensitivity to random density fluctuations, whereas a large j 
value will reduce that sensitivity at the cost of losing some 
locality. Through the use of Monte Carlo simulations Fer- 
reira & Lada (in preparation) find that a value of j — 20 is 
adequate to detect clusters with 10 to 1500 members. 

The positions of the cluster centres arc define d as the 
densi ty- weighted enhancement centres (jCasertano Sz HutI 
1 19851 ) 



where Xi is the position vector of the ith cluster member and 
Pj the jth NN density around this object. These centres do 
not necessarily correspond to the density peaks. 



3.2 Minimum spanning tree and Q 

The second method makes use of a minimum spanning 
tree (MST), a construct from graph theory, which is de- 
fined as the unique set of straight lines ("edges") connect- 
ing a given set of points without closed loops, such that the 
sum of the edge lengths is a minimum (|Kruskallll95^ : |Priml 
1957*). From the MST we d erive the mean edge length ^mst. 
Cartwri ght fc WhitworthI (l2004|) introduced the parameter 
Q ~ £mst/s, which combines the normalized correlation 
length s, i.e. the mean distance between all stars, and the 
normalized mean edge length ^mst- The Q parameter per- 
mits to quantify the structure of a cluster and to distinguish 



between clusters with a central density concentration and hi- 
erarchical clusters with possible fractal substructure. Large 
Q values (Q > 0.8) describe centrally condensed clusters 
having a volume density n(r) oc r~°', while small Q values 
(Q < 0.8) indicate clusters with fractal substructure. Q is 
correlated with the radial density exponent a for Q > 0.8 
and anticorrelated with the fractal dimension D for Q < 0.8. 
The dimensionless measure Q is independent of the number 
of objects and of the cluster area. The method has been 
successfully applied to bot h observed clusters and results 
of numerical simulations llCartwright fc WhitwOTthI |2004| : 
ISchmeia fclOes scn 200d; iKumar fc Sch mcia 2003). For de- 
tails of the method, in particular its i mplementation and 
normalization used for this study, see ISchmeia fc KlessenI 
(|2006t ). 



3.3 Application to the data 

The NN and MST statistical methods described above were 
applied to the c2d data of the Perseus, Serpens and Ophi- 
uchus molecular clouds. We first used the NN analysis (de- 
scribed in detail by Ferreira et al., in prep.) on the stellar 
and YSO populations in the entire clouds to detect den- 
sity enhancements and therefore embedded clusters. For this 
purpose, 20th NN density maps were generated for the full 
extent of the available c2d data. The 20th NN density is 
thought to be the optimum density to identify primary clus- 
ters and not pick up very loose groups (Ferreira fc Lada, 
in preparation). We found that using only the point sources 
that were classified as YSOs to produce the NN maps yields 
the best density contrast and boundary identifications. We 
determined the average 20th NN density in the clouds away 
from the clusters and define as cluster members all YSOs 
whose density exceeds 3a of the average value. The cluster 
boundaries are shown as thick lines in Fig. [TJ the density 
thresholds pthrcsh for the clusters are given in Table [2] 

All YSOs enclosed within the respective cluster defin- 
ing contour are treated as cluster members. For those we 
constructed the MST and determined Q. 

Note that the exact definition of the cluster boundaries 
is only relevant for the quantitative structure analysis using 
the Q parameter. Changing the cluster definition criteria 
may change the absolute Q values, but it does not affect the 
general trends and correlations discussed below. 



4 STRUCTURE ANALYSIS 
4.1 Cluster morphology 

The YSOs in the three investigated molecular clouds are 
found in the embedded clusters and in a distributed popula- 
tion that extends through the entire cloud (j0rgense n et al.l 
boodlHarvev et al.ll2007bl ). Figure [T] shows IRAC Channel 4 
(8 nm) images of the four largest embedded clusters IC 348, 
NGC 1333, Serpens cloud core and L1688, overlaid with the 
20th NN density contours. 

In the Perseus molecular cloud, apart from the em- 
bedded clusters NGC 1333 and IC 348, the region around 
Barnard 3 (B3) emerges as a moderately sized cluster, but 
no significant density enhancements are found associated 
with L1455 and L1448. Despite the large number of member 
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Figure 1. Spitzer IRAC Channel 4 (8 /im) images of the four main clusters IC 348, NGC 1333, Serpens A and L1688 overlaid with the 
20th nearest neighbour density contours. The thin contours are plotted in intervals of 40pc~^. The outermost thin contours correspond 
to a nearest neighbour density of 40 pc"'^. The thick contours indicate the cluster boundaries as defined in the text. The horizontal bars 
showing a projected length of 0.5 pc are based on assumed distances of 315, 260 and 135 pc to the Perseus, Serpens and Ophiuchus cloud, 
respectively. Crosses indicate the cluster centres, the known B stars are highlighted by diamond symbols. 



stars, the Perseus clusters show relatively small peak den- 
sities compared to the other clouds. Likewise, the Perseus 
cloud shows the lowest "background" YSO density of all re- 
gions. The density peak in IC 348 coincides with the B5 star 
HD 281159, the most massive star in this cluster, whereas 
the two B stars in NGC 1333 (BD +30 549 and SVS 3) 
are not associated with the densest region. The BO star 
HD 278942 is found close to the centre of the B3 cluster, 



however, it may be located behind the main Perseus cloud 
complex (|Ridge et al.ll2006l ). 

In the Serpens cloud two major clusters were detected: 
The main cluster in this cloud, usually calle d the Serpens 
cloud core, and referred to as Cluster A by iHarvev et al.l 
ll2006t). and a recent ly discovered cluster calle d Cluster B by 
iHarvev et all (|2006D and Serpens/G3-G6 by iDiupvik et all 
(200^). The Serpens core has a very steep density gradient 
and a prominent peak with a density of 1045 pc~^, the high- 
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Table 2. Clustering parameters of the identified clusters in the three regions. 



irtegion/ 




uec i^jzuuu^ 




Pthrosh 


Ppcak 






duster 


[deg] 


[ciegj 




fnc-^l 
[PC J 








Perseus 
















NGC 1333 


52.2563 


+31.3107 


2.09 


20 


261 


158 


0.81 


IC 348 


56.1151 


+32.1310 


2.53 


20 


140 


204 


0.80 


Barnard 3 


54.9612 


+31.9227 


1.67 


20 


45 


51 


0.69 


Serpens 
















Serpens core (A) 


277.485 


+1.2274 


1.00 


70 


1045 


84 


0.84 


Serpens B 


277.256 


+0.5059 


0.77 


70 


378 


43 


0.77 


Ophiuchus 
















Ophiuchus north 


246.348 


-23.4660 


0.78 


95 


126 


34 


0.68 


Ophiuchus centre 


246.374 


-24.1197 


0.52 


95 


213 


29 


0.73 


L1688 


246.683 


-24.5112 


1.85 


95 


509 


337 


0.72 


L1689S 


248.045 


-24.9294 


1.41 


50 


240 


78 


0.78 



est of all investigated clusters and twice the value of the 
second densest cluster. 

In the Ophiuchus cloud one highly structured, major 
embedded cluster was detected roughly centred on L1688, 
several smaller density enhancements clearly trace four fila- 
mentary structures that seem to converge on this prominent 
cluster. These density enhancements are in good correlation 
with the dust continuum emission maps. Two of these den- 
sity enhancements are picked up as small clusters by the 
NN method, we designate them Ophiuchus north and cen- 
tre. The Ophiuchus north cluster is centred around a group 
of bright B stars (p Oph D, p Oph AB, p Oph cf|. Further- 
more we detect the kn own cluster associated with L1689S 
iBontemps et al.ll200ll ). which shows a single peak and no 
distinct substructure. 

The shapes of IC 348, NGC 1333, Serpens A and L1689S 
can be roughly approximated by ellipses, whereas the hierar- 
chical L1688 cluster shows an irregular, nevertheless rather 
circular, perimeter. It is intersting to note that the elliptical 
clusters appear to have two peaks roughly coinciding with 
the foci of ellipses. However, these peaks seem to consist of 
of YSO populations of different evolutionary stages. This is 
seen in NGC 1333, with the density peak in the south and 
two luminous, massive B stars, which illuminate the promi- 
nent reflection nebula, in the north. In Serpens we note a 
group of bright stars north of the density peak, marked by 
a dotted circle in Fig.[T] The second peak in IC 348, seen to 
the south-east, seems to correspond to a new star formatio n 
episode close to HH211 (|Tafalla. Kumar fc BachilledbOOel ). 

Table [2] hsts all the detected embedded clusters along 
with their position, size, peak density, the Q parameter and 
the number of YSOs n* . The size corresponds to the length 
of the major axis of an ellipse fitted to the cluster areas. Its 
only purpose is to provide the reader with a rough estimate 
of the diameter of the cluster. The linear sizes are computed 
assuming distances of 315, 260 and 135 pc to the Perseus, 
Serpens and Ophiuchus cloud, respectively. 



^ These designations are not to be confused with the ones used 
for the subclusters discussed below and marked in Fig. [T] 



4.2 Q values 

Extensive analysis of the distribution of the MST edge 
lengths as a function of YSO class and the Q parameter 
as a function of YSO brightness was carried out only for 
four embedded clusters where the number statistics are sig- 
nificant and suitable for this kind of analysis. These four ma- 
jor embedded clusters are NGC 1333, IC 348, Serpens core 
and L1688. It can be seen from Figure [1] that the Serpens 
core, NGC 1333 and IC 348 clusters are centrally concen- 
trated, although at very different surface densities, whereas 
L1688 shows at least four peaks enclosed within the cluster 
defining contour. This is reflected in the Q values listed in 
Table [5] L1688 shows the lowest Q value (0.72) indicating 
fractal substructure, whereas the other clusters show Q val- 
ues of Q ^ 0.8, denoting central condensation. All the small 
clusters have values of Q < 0.8. 

In Table [3] we list the absolute and relative numbers 
of Class 0/1 and Class 2/3 objects and their Q values in 
the four large embedded clusters. Note that the classifica- 
tions (see Sect. [2} do not necessarily correspond to the stan- 
dard definitions. Assuming that the ratio of younger to older 
sources represents the relative ages of the embedded clus- 
ters, the data in Table [3] indicate that the L1688 cluster is 
the youngest of all, NGC 1333 and Serpens are slightly older 
and IC 348 is the oldest of all four. This is not identical, but 
quite similar to the age estimates based on more sophisti- 
cated methods mentioned in Sect. [T] In all the clusters but 
Serpens, the Q values are significantly higher for the Class 
2/3 objects than for the Class 0/1 protostars, indicating a 
more centrally condensed configuration for the older sources. 
In the centrally condensed Serpens A cluster the two groups 
show the same Q values. 

The three major subclusters in L1688 were investi- 
gated further, these are (following the nomenclature of 
iBontemps et al.ll200ll ') Oph A, Oph B and Oph EF (marked 
in Fig. [TJ. They are defined as having NN densities > 
180 pc~^, their clustering properties are listed in Tabled 
Oph A contains 41 Class 0/1 protostars and 26 Class 2/3 
objects and has a value of Q = 0.74 in the hierarchical range. 
Oph B (18 Class 0/1, 34 Class 2/3 sources) and Oph EF 
(16 Class 0/1, 27 Class 2/3 sources) show higher values of 
Q = 0.79 and Q — 0.77, respectively. Oph A is associated 
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serpens core 
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Figure 2. Normalized histograms of MST edge lengths for Class 0/1 (dotted line), Class 2/3 (dashed line) and all YSOs (solid line) in 
the four main clusters. 



Table 3. Numbers of objects and Q values for the different evo- 
lutionary classes 



Cluster 


Class 0/1 


Class 2/3 




n, Q 


n* Q 


NGC 1333 


62 (39%) 0.78 


96 (61%) 0.83 


IC 348 


72 (35%) 0.71 


132 (65%) 0.78 


Serpens A 


33 (38%) 0.85 


51 (62%) 0.85 


L1688 


194 (58%) 0.69 


143 (42%) 0.78 



Table 4. Clustering parameters of the subclusters in L1688 



Subclustcr 


RA 
[dcg] 


Dec 
[dcg] 


Ppcak 
[pc-2] 


n* 


Q 


Oph A 


246.553 


-24.416 


509 


67 


0.74 


Oph B 


246.836 


-24.465 


373 


52 


0.79 


Oph EF 


246.854 


-24.679 


367 


43 


0.77 



with a larger amount of gas and has a significantly higher 
Class 0/1 to Class 2/3 ratio than the other two subclus- 
ters, making it most likely the youngest of the subclusters, 
it shows the lowest Q value of the subclusters. Again we 
see an evolutionary sequence from the youngest, hierarchi- 
cal cluster to the more evolved clusters closer to central con- 
densation. When looking at Class 0/1 and Class 2/3 sources 
separately, we again find significantly higher Q values for 
the more evolved objects: In Oph A the value rises from 
Q = 0.74 for Class 0/1 to Q = 0.90 for Class 2/3, in Oph B 
from 0.77 to 0.81 and in Oph EF from 0.73 to 0.79. 



4.3 Distribution of MST edge lengths 

The lengths of the MST edges can be seen as the separa- 
tion from one object to its nearest neighbour. Histograms 
of the MST edge lengths indicate characteristic separations 
for YSOs in each of the clusters. Figure [2] shows the normal- 
ized histograms of the MST edge lengths in the four main 
clusters in bins of 0.02 pc, using the assumed distances given 
above (Section |4]2}. The histograms show a distinct peak at 
0.02 pc (Serpens, L1688) and 0.04 pc (NGC 1333), respec- 
tively, and the distributions spread up to separations of 0.2 
to 0.3 pc. In IC 348, however, two peaks appear at a larger 



separation of 0.04 and 0.08 pc. Also shown in Fig. [2] are the 
same histograms for Class 0/1 and Class 2/3 sources sep- 
arately. In the clusters NGC 1333, Serpens and L1688 we 
do not see a significant difference between the separations 
of Class 0/1 and Class 2/3 objects. Objects of all evolu- 
tionary states are more or less similarly distributed in the 
cluster. In IC 348, however, it is only the Class 2/3 objects 
which roughly follow the distribution of the overall cluster 
members, whereas the Class 0/1 objects seem to be homo- 
geneously distributed. This is because the central cluster in 
IC 348 is evolved and has blown out a cavity, with the cur- 
rent star formation occurring in a shell around the cluster 
( Tafalla et al.ll2006l : iMuench et "al]l2007l ). 

The characteristic edge lengths found from Fig.[2]do re- 
flect the true separations and are not a result of poor spatial 
sampling because these values are almost an order of magni- 
tude higher than the tolerance used for matching IRAC and 
MIPS data (4 arcsec, corresponding to 0.003 - 0.006 pc). We 
also investigated the MSTs and their edge lengths for the 
stars in the cluster held, but due to their large spatial den- 
sity they would generally yield very short edges and a narrow 
peak of the distribution at very small projected separations. 



4.4 Q versus magnitude 

Assuming that the IR source magnitudes are representative 
of the source masses and evolutionary stages we investigate 
the variation of Q as a function of magnitude. Q was calcu- 
lated for sources brighter than a certain magnitude in steps 
of 0.5 mag as long as the sample contained at least 40 ob- 
jects. This calculation was made for YSOs and stars sepa- 
rately. The resulting Q versus Channel 4 magnitude curves 
for YSOs and stars are shown in Fig. |3^ and b respectively. 
(Similar plots were made using Channel 1, 2 and 3 magni- 
tudes as well; they look similar in all bands.) The results for 
different clusters are shown using different line styles. Fig.|3K 
shows that the YSO curves remain flat or rising with mag- 
nitude for the clusters NGC 1333 and L1688 while they de- 
crease with magnitude for IC 348 and, less pronounced, for 
Serpens. This implies that in NGC 1333 and L1688 YSOs 
of all luminosities display a roughly similar configuration, 
whereas in Serpens and IC 348 the low luminosity YSOs are 
less centrally concentrated than the more luminous YSOs. 
In contrast, the stellar population shows a declining Q 
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Figure 3. Q as a function of source magnitude for YSOs (left) and stars (rigiit) in Channel 4 for IC 348 (solid line), NGC 1333 (dotted 
line), Serpens A (dashed line) and L1688 (dash-dotted line). 



versus magnitude curve for all the clusters, indicating that 
the fainter stars are relatively homogeneously distributed in 
space and the brighter stars are centrally concentrated. Since 
this curve is computed for a region encompassing only the 
dense cluster, the chances of picking up fainter background 
sources through the large column of extinction is low. Due to 
the same reason, stars that are cluster members outnumber 
foreground stars for the computed area. Therefore the ob- 
served effect is representative of the true variations among 
the cluster members. Furthermore, a random distribution 
expected from foreground or background objects will pro- 
duce a Q value constant over the whole magnitude range at 
Q « 0.73 and cannot reproduce the declining curve seen in 
Fig. Hj. 

Assuming the age of IC 348 as 3 Myr, that of Serpens as 
2Myr and that of NGC 1333 and L1688 around 1 Myr (see 
Sect. the YSO curves in Fig. |3^ show that the brighter 
objects appear to evolve towards centrally condensed con- 
figurations and fainter objects towards homogeneous distri- 
butions with time. The transition of the distribution of the 
lowest-mass/substellar objects from a centrally condensed 
configuration to a more homogeneous distribution is thought 
to be an effect of dynamical evolution. A detailed discussion 
of this effect in IC 348 and the Orion Trapezium cluster and 
its impli cations for the formation of brown dwarfs can be 
found in iKumar fc Schmeial l|2007l ). 



5 DISCUSSION 

5.1 Dynamical evolution and mass segregation 

The structure analysis of YSOs and stars presented in the 
previous section has several implications. The significantly 
higher Q values for the Class 2/3 objects compared to 
the Class 0/1 protostars show that the older sources a 
more centrally condensed than the younger objects, which 
in most cases show a hierarchical configuration. Also in 
other embedded clusters the Class 0/1 protostars are fo und 
in subclusters, e.g. in NGC 2 264 (iTeixeira et all liooel ) or 
the Orion Trapezium cluster l|Lada et al. I I2004I ). This is in 



agree ment with the predict ions of nu merical simulations 
( Bonnell. Bate fc Vine 2003. .; iSchmeia fc Klessen 2006) sug- 
gesting that clusters build up from several subclusters that 
eventually form a single centrally condensed cluster. How- 
ever, it is not possible to distinguish between the structure 
inherent to the cloud and that arising due to dynamical evo- 
lution. It is tempting to assume that the Class 0/1 objects 
trace the structure inherent to the cloud while the older 
objects represent the structure as a result from dynamical 
evolution. More than half of the mass of a Class protostar 
is found in its envelope. So a Class source can be seen as 
still part of the cloud, while the older YSOs are star-like ob- 
jects with a disc rather than an envelope and are therefore 
susceptive to dynamical interaction. 

While Q rises as a cluster evolves, its absolute value 
does not depend on the age of the cluster. NGC 1333 and 
IC 348 with estimated ages of 1 and 3 Myr, respectively, 
have similar Q values of 0.81 and 0.80, whereas the two ex- 
treme Q values in our sample (0.84 and 0.72) are associated 
with the clusters L1688 and Serpens thought to have ages of 
about 1 and 2 Myr, respectively. This shows that the cluster 
structure in general is not directly correlated with the evolu- 
tionary stage and rather depends on the existing large-scale 
structure of the cloud, but the observation that Class 2/3 ob- 
jects (which had more time and are also more susceptive to 
dynamical interactions) are more centrally condensed than 
the younger protostars is a clue that the cluster evolves from 
a hierarchical to a more centrally condensed configuration 
as predicted by numerical simulations. 

We interpret the decrease of Q with fainter (i.e. less 
massive) sources as an effect of mass segregation. However, 
the relation of luminosity to mass may be different for the 
different evolutionary classes, which may cause some inter- 
ference of evolutionary and mass segregation effects. The 
trend shown in Fig. [3^, based on Channel 4 (8 fim) mag- 
nitudes, is mainly representative of Class 2/3 objects, but 
it is similar in the other IRAC bands, and also in MIPS 
Channel 1 (24 fim), which includes Class 0/1 objects. 

A decrease of Q with fainter magnitudes is seen in 
IC 348 (~ 3 Myr) and, to a somewhat lesser extent, in 
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Table 5. Gas properties and Mach numbers 



Cluster 


Av 


T 




M 


Q 


Ref._^ 




[km s-i] 


[K] 


[km s"-'-] 








Serpens A 


0.77 


12 


0.21 


3.8 


0.84 


1 


NGC 1333 


0.95 


13 


0.21 


4.5 


0.81 


1 


L1688 


1.33 


15 


0.23 


5.8 


0.72 


2 



° Reference s for Av and T: 1; Ijiiina et al] lll999l') . 2: 
iMvers et al.] lll978h 

Serpens (~ 2 Myr) , whereas there are no signs of mass 
segregation in the YSO population of the youngest clus- 
ters. The stellar (i.e. older) population, on the other hand, 
shows significant mass segregation in all clusters. Therefore 
mass segregation seems to be a function of a ge. According 
to sim ulations of turbulent molecular clouds ( Bonnell et al.l 
I2OO3I ) YSOs form at different locations in the cloud and 
condense into a single cluster within a few free-fall times 
or less than 1 Myr. Dynamical interactions become impor- 
tant only after this initial assembly of YSOs is made. In 
particular, dynamical mass segregation occurs on a time- 
scale of the order of t he relaxation time of the cluster 
11 Bonnell fc DaviesHTooi ). This is consistent with the ob- 
served effects. Our findings are not necessarily contradict- 
ing claims that mass segregation is primordial, i.e. that 
massive stars form ne ar or in the centre of a cluster 
11 Bonnell fc Davieslll998h . In the youngest clusters, contrary 
to the older ones, the massive YSOs show a more hierar- 
chical distribution (smaller Q, see Fig. |3K) than the lower- 
mass objects, suggesting that the massive stars are still on 
the bottom of the potential wells of the former subclusters. 
This is consistent with theoretical scenarios in which merg- 
ers of small clumps that are initially mass segregated (or 
in which mass segregation is generated quickly by two-body 
relaxat ion) lead to larger clusters that are also mas s segre- 
gated (jMcMillan. Vesperini fc Porteeies Zwartll2007l ). 

The NN maps (Fig. ^ indicate that IC 348, NGC 1333 
and Serpens A have formed from a single dense core, while 
L1688 appears to lie at the interface of merging filaments 
(Ferreira et al., in preparation). L1688 has a highly hierar- 
chical structure and consists of at least three readily iden- 
tifiable subclusters, nevertheless it is detected as a single 
cluster by the 20th NN maps, whose centre does not coin- 
cide with any of the subclusters. Provided this geometrical 
centre corresponds approximately to the centre of mass and 
provided the density is high enough, the subclusters may 
merge to form a single centrally condensed cluster in the fu- 
ture. The structure of the cluster seems to be correlated to 
the cloud structure in the sense that the collapse of a single 
dense core more easily leads to a centrally condensed cluster, 
while interacting filaments are reflected in a longer-lasting 
hierarchical structure of the cluster. 



5.2 Turbulence and cluster structure 

As turbulence seems to play an important r ole in 



the build-up of a star cluster (lElmegreen et al.l 2000l : 
iMac Low fc Klessen|[200i : iBaUesteros-Paredes et al.ll2007h . 

we investigate the relation of the cluster structure with 



the turbulent environment. In Table [5] we list the veloc- 
ity dispersion Av and temperature T of the dense gas 
and the derived Mach numbers in Serpe ns, NGC 1333 and 
L1688 as traced by NH3 hue observations l|Mvers et al.lll978l : 
Ijiiina. Mvers fc Adams! [19991. The NII3 line is a tracer of 
dense gas which is rep resentative of the cores from which 
star clusters are born l|Bergin fc Tafallal |2007| ) and is also 
the tracer that it least affected by outflows. Therefore, the 
NH3 line widths better represent the turbulent properties of 
the cluster forming cores than other tracers of the molec- 
ular cloud such as CO or CS, which may be seriously af- 
fected by the presence of outflows and effects of mo lecular 
depletion (|Tafalla et al.ll2002l : [Bergin fc Tafallal 120071 '). Sim- 
ilarly, the best measure of the core temperatures is also 
obtained from the same emission lines. Hence we use the 
observed values of Av and T from this tracer to compute 
the Mach number M = Av/cs where the isothermal sound 
speed Cs = {RT/fi)^^^ with the specific gas constant R and 
the mean molecular weight jj, = 2.33. I Mvers et al.l (|l99ll ) 
argue that core linewidths are dominated by non-thermal 
motions (turbulence) only for core masses above ~ 22 Mq 
and that linewidths are thermal below ~ 7Mq. Since the 
core masses in all thre e clouds are abov e ^ 30 M© at an 
Av ~ 20 mag contour (|Enoch et al.|[2007i ). the assumption 
that the linewi dths are caused by turbulence is justified. 
The findings of iMvers et al.l (|l991 ) also lead to the conclu- 
sion that cluster-forming cores are turbulent w hile quiescent 
cores form isolated stars (see also I M ver3 1 200 ll ) . 

Comparing the Mach numbers and Q values in Table [S] 
shows that Serpens, being the most centrally condensed clus- 
ter, is associated with the lowest Mach number (turbulent 
energy) and L1688 being the most hierarchically structured 
cluster is associated with the highest Mach number (see also 
Fig. [l|. Note that we compare Q and M only for the three 
clusters that have a significant population of YSOs, are em- 
bedded in dense gas and have similar ages of 1 - 2 Myr. For 
this reason we excluded IC 348 from this comparison, which 
is significantly older and largely depleted of cold gas. Since 
Q seems to depend on the temporal evolution, a comparison 
only makes sense for clusters of roughly the same age. 

These results are in accordance with the scenario of 
turbulent fragmentation and hierarchical formation of clus- 
ters. Numerical simulations demonstrate that high Mach 
number fi ows lead to a highly fr a gmented den s ity str ucture 
( Ballcst oros-Paredes et al.ll2006l ). Enoch et all |20o3) study 
the same clusters and obtain a contradictory result; how- 
ever, the Mach numbers they compute are based on velocity 
dispersions obtained from observations of the CO molecule, 
which rather traces the outer envelope of the molecular cloud 
and may be contaminated by the effects of outflows. 

As this correlation is obviously based on small-number 
statistics and has not been found in numerical simulations 
of gravoturbulent fragmentation (Schmeja fc Kle ssen 200^ 
and considering the associated uncertainties of Av and T, 
it has to remain rather speculative. Nevertheless it is a hint 
that the structure of an embedded cluster may be somehow 
connected to the turbulent environment. Even though the 
turbulent velocity is probably not the only agent respon- 
sible for shaping the cluster, it may play a signiflcant role. 
While a high degree of turbulence in the cloud may keep the 
young stars in a more hierarchical distribution for a longer 
time, the absence of strong turbulent motions may help to 
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reach a centrally condensed configuration more quickly. As- 
suming that strong turbulence in the cloud also leads to 
a high velocity dispersion among t he embedded stars, this 
is also supported by the results of iGoodwin fc Whitwort"hl 
l|2004 ). who performed N-body simulations of the dynamics 
of fractal star clusters in order to investigate the evolution 
of substructure in recently formed clusters. They show that 
if the velocity dispersion in the cluster is low, much of the 
substructure will be erased, however if the velocity disper- 
sion is high, significant levels of substructure can survive for 
several crossing times. 



6 SUMMARY 

We used the Spitzer c2d survey data of the nearby molecular 
clouds Perseus, Serpens and Ophiuchus to identify embed- 
ded clusters and analyse their structures using the nearest 
neighbour and MST method. Apart from the large known 
embedded clusters IC 348, NGC 1333, Serpens core and 
L1688 we only found a few relatively small clusters. The 
Q parameter was determined for all embedded clusters in 
these clouds and the MST edge lengths and Q values were 
analysed as a function of YSO class and magnitudes. Our 
main results can be summarised as follows: 

(i) Among the four large clusters, IC 348, NGC 1333 and 
Serpens are centrally condensed (Q ^ 0.8), while L1688 
shows a hierarchical structure with several density peaks. 
The three main subclusters in L1688 show again a hierar- 
chical structure, with Oph A, presumably the youngest one, 
having the lowest Q value. The peak densities vary strongly 
between the clusters in different clouds, as does the density 
of the dispersed YSO population. 

(ii) In all clusters the Q values for the younger Class 0/1 
objects are substantially lower than for the more evolved 
Class 2/3 objects. This indicates that embedded clusters are 
assembled from an initial hierarchical configuration arising 
from the filamentary parental cloud and eventually end up 
as a single centrally condensed cluster, as predicted by nu- 
merical simulations of turbulent fragmentation. While the 
youngest, deeply embedded objects may represent the struc- 
ture inherent to the cloud, the distribution of the older ob- 
jects probably is a result from dynamical interactions. 

(iii) We find no signs of mass segregation for the YSOs in 
the youngest clusters, only for the YSO population in the 
presumably older clusters Serpens and IC 348. In contrast, 
the stellar population displays a clear mass segregation in 
all clusters. This suggests that the effect of dynamical inter- 
action becomes visible at a cluster age between about 2 and 
3Myr. 

(iv) The structure of a cluster may be related to the 
turbulent velocity in the cloud in a way that clusters in 
regions with low Mach numbers reach a centrally con- 
densed configuration much faster than those in highly 
turbulent clouds, in agreeme nt with N-body simulations 
l|Goodwin fc Whitwort"hll2004l '). 



The results from the Perseus, Serpens and Ophiuchus 
star-forming regions are in good agreement with the predic- 
tions of theoretical scenarios claiming that embedded clus- 
ters form from gravoturbulent fragmentation of molecular 



clouds in a hierarchical process. Our results may not be ap- 
plicable to the YSO populations in other regions, in par- 
ticular the Taurus star-forming region, which shows a level 
of clustering and velocity dispersions very different from the 
regions studied in this work and which is ge nerally harder to 
recon cile with the gravoturbulent scenario (|Froebrich et al.l 
l2006l) . 
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